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Abstract. We develop a scattering theory for perturbations of powers 
of the Laplacian on asymptotically Euclidean manifolds. The (absolute) 
scattering matrix is shown to be a Fourier integral operator associated 
to the geodesic flow at time it on the boundary. Furthermore, it is 
shown that on R™ the asymptotics of certain short-range perturbations 
of A fc can be recovered from the scattering matrix at a finite number of 
energies. 



1. INTRODUCTION 

In this paper, we develop a scattering theory for powers of the Laplacian 
on a class of manifolds which includes perturbations of Euclidean space and 
apply this theory to obtain new inverse scattering results on Euclidean space. 
This theory is a natural extension of the work of Melrose who devel- 
oped a theory of scattering for the Laplacian on asymptotically Euclidean 
manifolds. We show that the higher order scattering matrix has very simi- 
lar properties in this case. In particular, we show that the Melrose-Zworski 
calculus of Legendrian distributions, fll5l] , can be applied to construct the 
Poisson operator for the scattering problem and thus deduce that the scat- 
tering matrix is a Fourier integral operator associated to geodesic flow at 
time 7T. For higher order operators, this result appears to be new for the 
class of perturbations we consider even for R n . This theory is then also 
applied to extend the inverse results of Joshi and Sa Barreto on recovering 
the asymptotics of perturbations, ||, ||, [ 10 1 , (l^], to this higher order 
case. 

An asymptotically Euclidean manifold is a smooth manifold with bound- 
ary (X, dX) which is equipped with a scattering metric. A scattering metric 
is a smooth Riemannian metric, g, on the interior of X which blows up in a 
prescribed way at the boundary: there exists a product decomposition close 
to the boundary, p i— > (x, y) G [0, e) x dX, such that g takes the form 

dx 2 h(x,y,dy) 

g = — + 2 — ' 

x* x z 
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with h smooth on the closed space and hi x=0 a non-degenerate metric on 
dX. This is slightly different from Melrose's definition in E3] but was shown 



to be equivalent in The co-tensor, h\ x=0 , is independent of the decom- 
position chosen and thus we have a natural metric on dX. We shall assume 
throughout that a product decomposition close to the boundary has been 
chosen and fixed. 

It is important to realize that W 1 with the Euclidean metric is a special 
case of such a manifold. To see this, put x = \z\~ 1 and uj = z\z\ ; we then 
have 

o dx 2 dux 2 

which also shows that the induced metric on the sphere at infinity is the 
Euclidean metric on the unit sphere. 

Melrose showed that given a smooth function, /, on dX, and A G i \ {0} 
that there is a unique function u, smooth on the interior of X, of the form 

e iX/x x^f+ + e- iA/l x^/_, (1.3) 

with f± smooth functions on (X, dX) and /_ restricted to the boundary 
equal to / such that (A — \ 2 )u = 0. The scattering matrix is then defined 
to be the map, 5(A), on C°°(dX), defined by 

S(X) : / -> f +]dx . (1.4) 



It was shown in 14] that S(X) extends to a unitary operator on L 2 (dX) 
with the density induced by h. Melrose and Zworski ([^]) studied the micro- 
local structure of this operator and showed that it is a zeroth order, classical 
Fourier integral operator associated to geodesic flow at time 7r. In the special 
case that h(x, y, dy) = h(0, y, dy) + 0(x°°), Christiansen ([|]]) and Parnovski 



16|) showed that for A < 0, modulo smoothing, 



S(X) = ie"V^ x ^^. (1.5) 

Here we develop analogous results for perturbations of powers of the 
Laplacian. For k G N, we define a short range perturbation of A fc to be 
a symmetric differential operator of order 2k — 1 which close to the bound- 
ary can be written in local coordinates in the form x 2 P(x,y,x 2 D x ,xD y ) 
where P(x, y, r, rf) is smooth in (x, y) and is a polynomial of order 2k — 1 
in (r, rj). In particular, a real- valued, smooth function vanishing to second 
order at dX defines a short range perturbation. In Section § only, we allow 
certain pseudodifferential perturbations as well. 

We prove 

Theorem 1.1. Let V be a short range perturbation of A k and let A G M \ 
{0}. Then, given f G C°°(dX), there exists a smooth function, u, on X° 
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such that (A k + V — X 2k )u = and u is of the form 

e^x^U + e- iA /*x^/_, (1.6) 

with f± smooth functions on X such that /_ restricted to the boundary is 
equal to f. The function u is unique modulo smooth functions vanishing to 
infinite order at dX. 



The non-uniqueness here corresponds to the possibility of embedded dis- 
crete spectrum. The scattering matrix S(A) can then be defined precisely as 
before as the indeterminacy will not affect the lead term at the boundary, 
and we show that it has a unitary extension. 

We also prove 

Theorem 1.2. Let V be a short range perturbation of A k and let A E 
R \ {0}. Then 5(A) is a zeroth order classical Fourier integral operator 
associated to geodesic flow at time ir on dX. 



It is interesting to compare the scattering matrix for different powers and 
the following is an immediate corollary to our construction. 

Corollary 1.1. Let < k\ < k<i € N, and suppose Vj is a short range 
perturbation of A k i , of the form x l pj(x,y,x 2 D x ,xD y ), with pj(x : y : ^,r]) a 
polynomial in (£, f?). Let Sj(X) be the scattering matrix associated to A fcj + 
Vj — A 2fcj . Then Si (A) — 62(A) is a Fourier integral operator of order 1 — 1. 



It is therefore immediate that (1.5) holds also for the higher order scat- 
tering matrix when V = 0. 

In the special case of W 1 we study the problem of recovering asymptotics 
of a perturbation from scattering data. As in |12|, we need an aradiality 
condition to recover the perturbation. In fact, in [j^] it was observed that 
recovery is not possible without it. We shall say a perturbation is aradial 
modulo Schwartz functions if it is asymptotically equal to a sum 

00 

(1.7) 

1=2 a 

with each term fa,-iD" of the form \z\~ l times a composition of vector 

a 

fields tangent to the sphere. The aradiality conditions allows any zeroth 
order perturbation. 

Theorem 1.3. Let Vi,V2 be short range perturbations of A k on (M. n ,dz 2 ) 
such that V1 — V2 is an aradial differential operator of order I. Let Sj(X) be 
the scattering matrix associated to A k + Vj — X 2k , and suppose that for I + 1 
values of X > 0, Si (A) — S2(A) is smoothing. Then the coefficients ofV\ and 
V2 agree modulo Schwartz functions. 
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Note that we assume neither that V\ , V2 are of order I nor that they are 
aradial. 

Our approach to this higher order scattering problem is highly influenced 
by that of Melrose, [O], and that of Melrose- Zwor ski, p5| . In particular, 
we use the scattering calculus developed by Melrose and used for the case 
k = 1 to study the general case and use techniques similar to those of [14] 



to establish the existence of the scattering matrix. To establish the micro- 
local structure of the scattering matrix we proceed as in [15| to construct 



the Poisson operator for the scattering problem as a Legendrian distribu- 
tion associated to a pair of intersecting Legendrian submanifolds using the 
calculus developed there. 

To prove the inverse result, we choose to follow the approach of || rather 
than [^] in order to increase the readability of the paper for non-experts. 
In particular, we establish our results for the case of R™ without explic- 
itly using the Melrose- Zwor ski Legendrian calculus. The reader may regard 
these proofs as a warm-up for the construction of the Poisson operator in 



the general case. As in ||, |J, [11], [12], the proof proceeds by estab- 



lishing that the principal symbol of the difference of the scattering matrices 
determines and is determined by a weighted integral of the lead term of the 
difference of the perturbations over geodesies of length tt. The injectivity of 
this transformation is then deduced by using some elementary calculus and 
some deep results of Bailey and Eastwood, ||, on the integral geometry of 
tensor fields on projective space. 

The main reference for higher order scattering on W 1 is Chapter 14 of |?| 
where a scattering theory for perturbations of a much more general class of 
constant coefficient operators on W 1 is developed - this extends ideas devel- 
oped by Agmon and Hormander in Q and Agmon in Q. The complementary 
problem of studying the recovery of compactly supported perturbations of 
higher order operators on W 1 was studied by Liu in | |13| |. 

The results of Bailey and Eastwood ( || ) which we use in section ^| to prove 
the inverse results are restricted to projective spaces and spheres. Since in 
addition, we believe the principal interest in this problem is in Euclidean 
space, we restrict ourselves to studying that important special CcLSG ELS the 
analysis is much more accessible. The problem of recovering the asymptotics 



of metrics for the k = 1 case has been studied in [11] and it is likely that a 



similar result could be proven in the higher order case. 

The question of recovering the entire perturbation from the scattering 
matrix at fixed energy is an interesting one but is still open even in the case 
with k = 1. 

We are grateful to Mike Eastwood for helpful conversations and to the 
Royal Society for the travel grant which made those conversations possible. 
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2. Basic scattering theory of A k + V 

In this section, we construct solutions of A k + V — X 2k having specified 
behaviour at the boundary, leading to the definition of the scattering matrix. 

Throughout, A denotes the Laplacian associated to a scattering metric 
on a compact manifold (X,dX), k is a positive integer, and A G R \ {0}. 
We assume a product decomposition with boundary defining function x has 
been fixed. We call the smooth functions vanishing to infinite order at the 
boundary the Schwartz functions on X. We use the product decomposition 
to extend functions on the boundary smoothly into X by making them 
constant in the normal direction and cutting off. 



The results of this section are very closely related to the results of [14]. 
We have tried to state the results in a manner accessible to those unfamiliar 
with that paper, but in order to avoid repetition we omit proofs which 
follow essentially as in |l4j , only giving some indication of how they need to 
be modified to work in this setting. We recall some of the definitions and 



results of [14 1 but refer the reader to the original for full details. 



In this section only, we will allow a somewhat larger class of perturbations 
of A k . In order to describe them, we recall some notation. 

Let V SC (X) be the space of all smooth vector fields of finite length with 
respect to a scattering metric. Near a point on the boundary, x 2 -^ and 
a 'c5y 7 ' * = n ~ 1> f° rm a basis for V SC (X), where j/, are coordinates on 
dX. The set of scattering differential operators of order m is 

Diff™(X) = span < J < m (V sc (^))l. 

The short range perturbations which we are allowing in the remainder of 
this paper are elements of x 2 Diff^^X). 

Most of the results of this section hold for a wider class of perturbations: 
elements of ^ 2k ~ 1,2 (X), part of the (small) calculus of scattering pseudo- 
differential operators. We refer the reader to Jl4] , Sections 4 and 5] for the 
full definition and some properties of (X). However, we remark that 
^Diff^(X) C 9^(X) and (A - z)' 1 G *- 2,0 (X) when z # [0, oo) (@, 
Theorem 1]). Moreover, if we consider the manifold (M. n ,dz) (which by an 
earlier discussion under radial compactification becomes a compact manifold 
(M n ) rc with scattering metric), the operators corresponding to 
have Schwartz kernels of the form 



A(z,z') = j^ J e^<a L (z,()d( 
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with cll satisfying 

\D%D$a L (z,0\ < C a p(l + + \Q\) m ~ m 

(0, (4.1),(4.2)]). 

Definition 2.1. We shall say a differential operator V is a short range 
perturbation of A k if it is symmetric and ifV€ x 2 Diff 2 ,? _1 (X). We call a 
pseudodifferential operator a pseudodifferential short range perturbation of 
A k if it is symmetric and it is an element of ty 2k ~ 1,2 {X). 

Throughout this section V is a pseudodifferential short range perturbation 
ofA fc . 

We prove 

Proposition 2.1. Let V be a pseudodifferential short range perturbation 
of A k and A G R\ {0}. Given f G C°°(dX), there exists u G C°°(X), 
unique modulo Schwartz functions, such that u\qx = / and (A k + V — 

\ 2k )(e~ lX / x u) is Schwartz. 

Note we make no assumptions on the sign of A. 

This will follow easily once we have proven 

Lemma 2.1. Let V be a pseudodifferential short range perturbation of A k 
and AeR\{0}. // / eC°°(dX), then 

(A k + V- \ 2k ){e- lX /' x x^ +a f) = e~ iX / x {k\ 2k ~ l C a x^ +a f + x^ +a g) 
with g G C°°(X), where C = and C a ^ for a ^ 0. 



Proof. For k = 1, this follows as in using [15, Lemma 8] for the mapping 



properties of pseudodifferential short range perturbations. For general k, one 
simply iterates. □ 



Now to prove Proposition |2.1| , we simply choose the lead term to be / 
and then repeatedly iterate to compute all the terms of the Taylor series. 
The result follows from Borel's lemma. The uniqueness follows from the fact 
that C a / for a + 0. 

We will repeatedly use the operator 

Q = Q(A) = ^ A^A^'- 1 (2.1) 
3=0 

since A fc — \ 2k = Q(A — A 2 ) = (A — X 2 )Q. Using the same techniques as 
Theorem 1 of 0, one can show that for A G M\{0}, Q(A)" 1 G * sc 2fc+2 '°(X). 

A key idea of many of the proofs is that 

A k -X 2k + V = Q{A-\ 2 + V), (2.2) 
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where V = Q V £ ^l^(X). The importance of this is that the symbol 
of V vanishes to second order at the boundary. As noted in Remark 3 



of [14], many of the results there hold if A is replaced by A + W, when 
W £ x 2 C°°(X), because they depend on the properties of the principal 
symbol and the boundary symbol, which are unchanged by the addition of 
such a W. What we are doing here is allowing a somewhat more general 
perturbation, but with the same kind of decay at the boundary. Thus, the 
results of Melrose's Propositions 9-11 hold if A is replaced by A + W, with 

The following proposition is closely related to Proposition 11 of |14| and 
follows from the modifications of Propositions 9-11 indicated above and from 

Proposition 2.2. If X € M\{0}, u £ L 2 SC (X), V is a pseudodifferential short 
range perturbation of A k and (A k + V — X 2k )u = 0, then u is Schwartz. 

If k = 1 and V is a short range (differential) perturbation, then A + V has 
no positive eigenvalues. However, for (differential) short range perturbations 
of A fc , k > 1, the situation can be different. 

Proposition 2.3. There are (differential) short-range perturbations V of 
A k , k > 1, even, such that A k + V has positive eigenvalues. 

Proof. We give two ways of constructing such examples. Let j € N. 

Using the min-max principle, for any asymptotically Euclidean manifold 
X one can find V € x 2 C°°(X) so that A + V has a negative eigenvalue —a. 
Then (A + V) 2 ^ has a positive eigenvalue a 2 K It is possible to choose V to be 
compactly supported and such that A + V has many negative eigenvalues. 

To construct a Schwartz potential V such that A 2j + V has an embedded 
eigenvalue, choose r > and let x be a globally defined boundary-defining 
function, x € C°°(X). Then 

(A 2j - T 4j )e~ T/x x (n ~ 1)/2 = e>(§(\+ 9 )/ e ]- T /§). 



Just as in Proposition |2,1| we can use Lemma |2.1| to successively solve away 
the errors, resulting in a u which satisfies 

(A 2 J - T 4j )u = g = 0(§°°l- r/§ ) 

with u = e~ T / x x( n ~ 1 ^ 2 (l + C(§)). In fact, by doing the asymptotic sum- 
mation judiciously, we can ensure that x^^ n ~ l ^ 2 e T / x u is nowhere vanishing. 
Set V = —u~ 1 g. Then A 2 - 7 -\-V has eigenvalue r 4 - 7 with eigenfunction u. □ 

We continue with the results needed to define the scattering matrix. 

Theorem 2.1. If X £ R \ {0}, V G ^ 2 ( ?~ 1 ' 2 (X) is symmetric, and f is a 
Schwartz function on X , orthogonal to the L 2 null space of A k + V — X 2k if 
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there is any, then there exists u with x ( n l >' 2 e l>i ' x u G C°°(X) such that 

(A k + V- \ 2k )u = f. 

This is an analogue of Propositions 12 and 14]. 

Proof. As in Proposition 14 of |l 



u = lim(A fe + V - (A + itf k r l f 

no 

where the limit exists in x -1 / 2-5 H^(X) for any 8 > 0. To show that the 
limit exists, let u t = (A k + V - (A + it) 2k )~ l f. Then, since / is orthogonal 
to the null space of A k + V — X 2k , so is Uf for i > 0. Having made this 
observation, the proof follows as in Proposition 14], giving us the same 
additional microlocal regularity as well. 



To finish the proof, one uses the analog of Proposition 12 of Hj. That is, 
if u G C-°°(X), WF*r 1/2 (u)nR + (X) = 0, and (A+Q^V-X^u G 
then e -*A/x x -(n-i)/2 { 2 G C°°(X). The proof follows just as the proof of @, 
Proposition 12], first noting that [Ao,^ 1 ^] G ^s^PO and thus 



The second observation that allows the proof to proceed as in [14, Prop. 
12] isthatifDiff£(A> C H™> m (X), then Diff^'(X)ADiff^(X)u C H^' m+s (X) 
when I' < I is a nonnegative integer and A G ^^f(X). This then allows the 
inductive step in the proof of p3, Proposition 12] to proceed just as it does 
there. □ 

We shall also use the following boundary pairing result, the analogue of 



[14, Proposition 13]. 



Proposition 2.4. Suppose V G ^ 2k 1)2 {X) is symmetric, X G R \ {0}, 
Ui, i = 1, 2 satisfies (A k + V — X 2k )u{ = fi, fa is Schwartz, and = 
x (n-l)/2( e i\/* a + + e -*V* a -) ) W uh af G C°°(X). Lei &f = {af)\ dx - Then 

2ikX 2k - 1 ( (bfbf - bfb^)dh = ( u{fe - fmdg. 
Jdx Jx 

Proof. Using Q = Q(X) = J2jZo X 2j A k ~ j ~ l , we have 

/ uif 2 - fiu^dg 
Jx 

= I ui(A - A 2 + VQ~ 1 )Qu 2 - Q(A - A 2 + Q- x V)u x uidg 
Jx 

= I ui(A-X 2 + VQ- 1 )Qu 2 - (A - A 2 + Q-^mQu^dg 

Jx (2.3) 
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as Q is self-adjoint and (A - A 2 + Q~ l V)ui G C°°(X). Recall that V is 
self-adjoint as well. Choose a function (j> G C°°(M) such that <f>(t) = if 
t < 1 and = 2 if t > 2. Then (2.3) is equal to 

lim / 0(x/e) («i(A - A 2 + VQ~ 1 )Qu 2 - (A - A 2 + Q -1 ^)«i^) d» 
40 Jx v y 

= lim / [A + Q-V, ^(x/e)]in^d^. 
40 Jx 

Since V is short-range, lim e iQ J x [Q~ 1 V,4>(x/e)]uiQu2dg = 0. To compute 
the remainder, we just use the corresponding results from [14, Proposition 
13] and the fact that at the boundary of X, 

□ 



We will need 

Theorem 2.2. //Asl\ {0} ; V is a pseudodifferential short range pertur- 
bation of A k , and f G C°°(dX), then there exists f± G C°°(X) such that 
u = x^ (e lX l x f + + e~ iX / x f-^ satisfies 

(A k + V- \ 2k )u = 

and the restriction of /_ to dX is f. The function u is unique if X 2k is not 
an eigenvalue of A k + V and is unique up to the addition of a Schwartz 
function if X 2k is an eigenvalue. 



Proof. Most of the proof follows from Proposition 2A_ and Theorem 2.1. We 
use Theorem 2J. to solve away the error obtained by constructing the formal 
expansion as in Proposition |2.1| . We note that since the error we wish to 
solve away is of the type (A k + V — X 2 )g, where g G x~ 1 ^ 2 ~ e L 2 c (X), and since 
eigenfunctions are Schwartz we may still integrate by parts to obtain that 
the error is orthogonal to the eigenspace. Finally, suppose that there are 
two such u satisfying the conditions of the theorem. Then their difference 



V — X 



2k\ 



and v 



„(n-l)/2 iX/x 



g + + x 



(n+l)/2 p -iX/x 



9- 



v satisfies (A k 

with g± G C°°(X). Then the boundary pairing result of Proposition 2.4 
gives us that (g+)\dx = 0, and thus v is in L 2 C (X). Consequently, v is an 
eigenfunction and is thus Schwartz. □ 



We can therefore define the scattering matrix as in the usual case: 

Definition 2.2. The scattering matrix is a map on C°°(dX) taking a func- 
tion f to /+ restricted to dX where /, /+ are as in Theorem \2.2\. 



Using Proposition 2.4 and the uniqueness it implies, it is easy to check 
that the scattering matrix can be extended to a unitary operator on L 2 (dX) 
which satisfies S^A) -1 = S(— A). 
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We remark that if V € ^f 2 ^~ 1,2 (X) is symmetric, then A k + V, initially 
considered as an operator on the Schwartz functions, has a unique self- 
adjoint extension to a domain in L 2 C (X). We end this section by describing 
the continuous part of its spectral measure. Let P be the Poisson operator 



whose existence is given by Theorem 2.2 



P{\) : C°°(dX) 3 /h U £ x { - n - 1 ^ 2 e iX / x C co {X) + x^-^V^C^pT) 

where u is the function in Theorem [T^. Then the continuous part of the 
spectral measure of A k + V is given by 

dE c (X) = ^-P(X)P*(X)dX 

where X 2k is the spectral variable and A € [0, oo). This can be seen by first 
noting that if V = this follows from |[|, Lemma 2.2] (see also || Lemma 
5.2]) and the general case then follows as in [17, Appendix to XL 6]. We note 
that A k + V has no singular continuous spectrum, though, as previously 
noted, there may be discrete spectrum. 



3. The Poisson Operator for Euclidean space 

As an introduction to the general case, we use the techniques of || to 
construct the Poisson operator in the Euclidean case, where the construction 
can be done more explicitly. Following Q we look for a Poisson operator 
in the form e lXz ' UJ a(z, uo) with a a poly homogeneous symbol in z and u is a 
smooth parameter. 

Let a(z,u) be any polyhomogeneous symbol in z, with lo a parameter. 
Then 

A(e iA2XJ o) = (A 2 a - 2iAcu.^ + Aa)e iXzM , (3.1) 

so we conclude that 

A k (e iXz - UJ d) = (X 2k ti - 2kiX 2k ~ l uj.— + 6)e iAz - w , (3.2) 

ox 

with b a symbol two orders lower than a. 

Now V is a short range perturbation of A k , which for W 1 means that 

V= f« D *> ( 3 - 3 ) 

|a|<2fc-l 

with f a G S~ 2 (M. n ). We conclude that it maps e iXz - u S m to e iXz - u S m ~ 2 and 
hence that if a S S™ hg we have, 

(A k + V - A 2fc )(e iA2 - w a) = (-2HA 2fe - 1 cu.^ + c)e iXz ^ ', (3.4) 
withcGS-" 2 . 



HIGHER ORDER SCATTERING 



11 



Returning to the Poisson operator, by taking the lead term of a to be 
the constant function 1 we have an error in . To solve away this and 
subsequent error terms, we proceed precisely as in || Section 2], using a 
term in S~i g +1 to solve away an error term in S p £ g (modulo terms in S~£ g 1 ). 
We solve the transport equations, 

-2ki\ 2k ~ 1 u;.^ = d, 
oz 

where d is the error, along the geodesies on the unit sphere from u to —oj. The 
only difference from the k = 1 case of |^] is a factor of A 2 ^ 1 ). As before, the 
transport equation degenerates on approach to —to and in particular blows 
up like (tt — s)~ r as s, the geodesic distance from u>, tends to it, when solving 
for the term of homogeneity of order — r. 

So as before we need a different ansatz close to — uj. We can locally in u 
rotate our coordinate system so that w is the north pole. So close to the 
south pole, we look for the Poisson operator in the form 

S a e i\(Sz'.»-VWsS\z\) a (^-,S,^j dSd/X, (3.5) 

with a(t, S, /i) a smooth function compactly supported on [0, e) x [0, e) x S n ~ 2 . 
We denote the class of functions that can be written in this form plus a 
Schwartz error by J 7 ' Q . We recall from ||] that this class is asymptotically 
complete in 7. It follows from a stationary phase computation carried out in 
H that away from the south pole this is equivalent to the previous ansatz 
with a symbol of order — 7 + ^2 ■ 

We recall from Q, 

Proposition 3.1. Ifu(z,u) £ H' a and f E C°°(dX x dX) then 

e a l*l j u{\z\6,oj)f{e,uj)d6duj 

is a smooth symbolic function in\z\ of order — 1 — a and its lead coefficient is 
\z\~ a ~ l (K, f) where K is the pull-back of the Schwartz kernel of a pseudo- 
differential operator of order a — 7 — (n — 2) by the map 9 h —9. The 
principal symbol of K determines and is determined by the lead term of the 
symbol, a(t, S, ^), of u as S — > + . 




We also have as a special case of Proposition 15 from [15] 
Proposition 3.2. If u £ J 7 '*, then 

(A k - X 2k )u £ P +1 ' a+1 , 

and 

vu £ r+ 2 ' Q + 2 . 



It is also important to note that, 
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Lemma 3.1. If u € I 7 '^ and (A k + V - \ 2k )u e JT+^V, then (A k + 

7-A 2 *)«€iT ,+i ' 2 * 1 . 

This follows from a slight extension of the argument in the proof of Lemma 
3.2 in §. 

We also recall from Q that 
Proposition 3.3. If u € = fV 7 '", ^en u = e-* A l z l/(z), raitt / a 

7 

classical symbol of order —a — 1. 

To carry out our construction we first use the original ansatz, obtaining a 
symbol which blows up on approach to the south pole. Near the south pole, 
we use the second ansatz, ( |3.5|) , taking 7 = — 2=- ; a = 2=2. The argument is 
then identical to the one in |J. At the end of the construction, we obtain an 
approximate Poisson operator, P, such that (A k + V — X 2k )P(X) is Schwartz. 
This error will in fact be orthogonal to the eigenspace at energy A (if there 
is one) as the eigenfunctions are necessarily Schwartz and the orthogonality 
follows from self-adjointness and integration by parts. We can therefore 
apply the resolvent (Theorem 2.1) and gain a term of the from e jA ' 2 '/ with 



/ smooth, solving away the error. 

4. The inverse problem 



In this section we prove Theorem 1.3. Suppose we have two short range 



perturbations, V\, V2, of A fc and we want to compare their scattering matri- 
ces. In particular, suppose 

V 1 -V 2 = a *(z)D°, (4.1) 

[a|<2fc-l 

with a a G Sphg~ r Q^ n )- If we carry out our construction for each Vj, then 
the first r terms of the construction with the first ansatz will be the same 
but the term of homogeneity — r will be different. In particular, the forcing 
terms in the transport equations will differ by 

W- r = e- iXzM (Vi - V 2 )e lXz - UJ = J2 w a a a _ r _iAl a l, (4.2) 

|a|<2fc-l 

where the lead term of a a . We therefore conclude that the terms 

of homogeneity — r will differ at (7(5), uS) by 
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2/cA 2fe - 1 (sins) 7 



W- r (-f(s'),uj)(sms') r - 1 ds', (4.3) 







where 7 is a geodesic on the sphere running from u to —uj. As we have 
shown that the lead singularity of the first ansatz as s — > it— is essentially the 
principal symbol of the scattering matrix, we conclude by the same argument 
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that the difference of the scattering matrices will be of order — r and that 
the principal symbol of the difference will determine and be determined by, 



1 



A 2fc-1 



W_ r (7(s), u) (sin s'Y^ds, (4.4) 



o 



for all geodesies 7 with u equal to 7(0). 

Note that as W depends polynomially on A, so does the principal symbol. 
In particular, if we know the scattering matrix for 2k different values of 
A > 0, then we can separate out the parts coming from the differing orders 
of \a\ in the forcing term. More generally, if we have a perturbation of order 
I < 2k — 1, then we determine the asymptotics from I + 1 values. So in the 
sequel, assuming the perturbation is of order I < 2k — 1, we only consider 
forcing terms of the form 

E ^(z)D», (4.5) 

where I < 2k — 1. To solve the inverse problem, then, we must show that 
a a (z), \a\ = I, can be recovered from knowledge of the transform ( |4.4| ), where 
we replace the sums in ( |4.lD and (|4.2|) by sums only over the terms with \a\ = 
I. That is, we need to show that the transform ( |4.4j ), a weighted integral 
along geodesies of length tt, is invertible. The injectivity of the transform 



was proven in [10] for the case I = 0. For higher order perturbations (Z > 
0) the question is more subtle and requires the notion of aradiality. We 
say a perturbation is aradial if it has no radial component, that is, if the 
perturbation is a span of vector fields tangent to the sphere. Clearly, all 
zeroth order perturbations are aradial. The injectivity of the transform 



( ^4.4[) for aradial first order perturbations was shown in 12] and for second 



order aradial perturbations in [11]. The impossibility of recovering first 



order perturbations that are not aradial was also shown in [12]. We now 
look at the general case. 

First we re-express the forcing term more invariantly. We regard the 

perturbation as an Z-form, \x = a a - r -i{z)dz a . Since dz a is symmetric, 

\a\=l 

/i is symmetric. The aradiality means that /U makes sense as a form on the 
sphere and is determined by its values as a map from the tangent space of 
the sphere to R. 

To re-express the forcing term we rotate so that to = (0, . . . , 0, 1) and 
7(s) = (0, . . . , 0, sin(s), cos(s)). We then have -£;{s) = (0, . . . , 0, cos(s), — sin(s)). 

7T 

The transform is then just / a( 0i ... i o,i)(7('S))(sin s) r ds. We show that this 



is equal to 



(-ly^sin^- 1 ^,^^)^^ 
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where (fj,, ^(s)/ is the pairing of // with ^r® - "®^: copies). To prove this, 
observe that as 7 lies in a plane this is really a two-dimensional question, so 
without loss of generality its enough to take n = 2. Now fj, is a symmetric 
/—form so 

i 

3=0 

Thus 

0*. = E «(jW)(7W)(- sin( S ))^(cos( S )y . 

i=o 

It follows from aradiality and the fact that the perturbation is of order I 
that 

A / — cos(s)V 



and thus we have 

(- s in( S ))V^(s)> =o ( o,i)E ^Vcos( S ))^(sin( S )) 2 ^, 

j=o V/ 

which is of course equal to a(o,z) • This establishes the equality. 

We now want to show that if /i is a symmetric /—form on the sphere and 

(/ifrOO), ^( S ))(sin( S ))^- 1 d S = 0, (4.6) 


for all geodesies 7 on the sphere, then fi = 0. Since we know this is true for 
any geodesic 7, we have 

Il+r-l^a = J 0*(7(s + ^f( s + a))(sin(s)) /+r - 1 ds = 



for any a. Differentiating with respect to a (see [|10[), we deduce that Ij t <y, a = 
implies that Ij-2,-y,a = and thus, differentiating repeatedly, if / + r — 1 
is even and I/ +r _i i7jQ = then \i is even. Also, if / + r — 1 is odd and 
Il+r-i,7,a = 0; then fi is odd. 

Since we can always reduce by two, it is enough to consider the cases 
where r is 1 or 2. Considering r = 1, we have for any geodesic 7 that 

7T 





If we take a geodesic starting at z n = 0, then we deduce that 
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By rotational invariance we have that 

7T 

W)( 7 ( S )),^(s))d S = 0, 
o 

for all homogeneous polynomials, p, of order I. As p\x is even, we can regard 
it as a symmetric tensor on projective space which is in the kernel of the 
generalized x-ray transform on symmetric /—tensors. Fortunately, the kernel 
of this operator has been identified by Bailey and Eastwood, ||. They 
showed that the kernel is precisely the symmetrized covariant derivatives 
of symmetric (I — l)-tensors. Therefore, to show that [i = we need to 
show that if p[i is a symmetrized covariant derivative for all homogeneous 
polynomials p of order I, then [i = 0. 

We show this by showing that \x has to vanish at the north pole. By 
rotational invariance it will then follow that \x vanishes everywhere. For the 
case r = 2, we apply the same arguments to z n [i and the result will follow 
from the case r = 1. 

We first prove 

Lemma 4.1. Let \x he a symmetric co-tensor of order I on M n_1 such that 
pfj, is a symmetrized covariant derivative for every homogeneous polynomial 
of order I. Then fi vanishes at the origin. 

Proof. To prove this lemma, we introduce some new notation. If a = 
(ai, . . . , a r ) € {1, . . . , n — l} r then 

Note that d a ^ d a even when both sides make sense. We also put 

dx a = dx ai . . . dx ar . (4.8) 

Let ott = (ai, . . . , a t -i,a t+ i, ■ ■ ■ a t ). 

The symmetrized covariant derivative, V s r/, of a tensor r\ is obtained by 
taking the usual covariant derivative and then averaging over the symmetric 
group to make it symmetric. On M" --1 , if the symmetric I — 1 tensor r/ = 
^(fr^dxry, and V<j7/ = ^ i> a Ax a , then 

1 1 

3=1 

Our proof of the lemma follows from the observation that the symmetric 
/ — 1 tensors must satisfy certain PDEs. If a, /3 £ {1, ... ,n— l} r , we define 
an exchange to be a map swapping certain places in a with ones in (3. We 
include the case where no swaps take place. The order of the exchange is the 
number of swaps. The sign of the exchange will be —1 to the power of the 
order. If the exchange is e, we denote the sign of e by sgn(e) and the new 
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value of a by e(a,f3)^ and of (5 by 
cki with Pi then 

e(«,/?) (1) = 

e(«,/3) (2) = 
sgn(e) = 



e(a, P)^ ■ For example, if e exchanges 

(Pi,a 2 , ... ,a r ) 

(ai,p 2 ,...,P r ) 
-1. 



If a, P € {1, . . . , n — l} 1 , and \i = YH J 'idx 1 is a covariant symmetrized 
derivative then 

^sgn(e)a e(ai/3)(2) V e(ai/3) (i) =0, (4.10) 

e 

where the sum is taken over all exchanges. To prove this one substitutes the 
expression for a covariant derivative and observes that the terms with no 
swaps will cancel with some of the terms from swaps of order one but that 
the swaps of order one will then have remainder terms which are canceled 
by swaps of order two and so on. 

Suppose we have that every homogeneous polynomial of order I times \i 
is a symmetrized covariant derivative. We show that each term tp a vanishes 
at the origin. We consider two cases. The first case is that a does not 
contain all possible values from {1, . . . ,n — 1}, (which will always happen 
when I < n — 1). Let r be the value not taken. Let P = (r, . . . , r). We have 
that 



E 



s g n ( e )d e ( Ql/ 3)(2) (&/9^e(a,/J)W ) = ( 4 - n ) 



where xp = xr x xq 2 ■ ■ ■ xr v Upon evaluating at x = 0, all terms will vanish 
except da(x8ip a ), which will equal ifj a , which must therefore be zero. 

In the second case, a takes all values from {l,...,n — 1}. We let f3j 
equal 2 if ay equals 1 and 1 otherwise. As before, ( [4.11 ) is satisfied. Upon 



evaluation at zero, the only terms that will not vanish are those for which 
e(a,P) = (a, P) (up to the order of otj and the order of Pj) and such e will 
have positive sign as an even number of swaps will be involved. We therefore 
have that a positive multiple of ip a is zero and thus that ip a is zero. □ 



This completes the proof for M n with the Euclidean metric. We, how- 
ever, wish to obtain a similar result for the sphere. Projecting the upper 
hemi-sphere onto the plane, we obtain coordinates on the sphere and the 
metric agrees with the Euclidean one at the north pole (which corresponds 
to the origin). The symmetrized covariant derivative will then agree with 
the Euclidean one up to terms vanishing there. As we are applying the result 
to polynomials of order / these additional terms do not affect the argument 
and so we conclude that the result will also hold for the sphere. 
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So we have shown that if the principal symbol vanishes then the lead term 
of an aradial perturbation vanishes also and thus by induction the asymp- 
totics of aradial perturbations are recoverable from fixed energy scattering 
data. 



5. Review of Legendrian Distributions 

In this section we review and rephrase the material we need from [14| and 
[15|. Here X is a compact manifold with boundary dX and g is a scattering 
metric on X and we have chosen a product decomposition of the for m (|1.1| ) . 
Our account is necessarily brief and we refer the reader to [14] and [p.5[| for 
more details. 

There is a natural bundle over X called the scattering cotangent bundle 
which is denoted SC T*{X). This is the dual to the bundle of smooth vector 
fields of bounded length with respect to some (and hence all) scattering 
metrics on X. The restriction of SC T*(X) to dX is denoted sc T*{X)\ dx and 
carries a natural contact structure. If y are local coordinates on dX and /x 
are the corresponding dual coordinates, then (y, //, r) form local coordinates 
on sc T*{X)\ dx , where r is the coefficient of 

We recall from Section || that a differential operator P(x,y,xD y ,x 2 D x ) 
will be in SC Q 1 / 2 ) if it is of order m and the total symbol as an 

operator in xD y ,x 2 D x vanishes to k th order at the boundary. It then has a 
well-defined symbol at the boundary 

j(P) = x k Pk + x k+1 Pk+1 , Pk ,p k+1 E C°°(R x T*8X). 

Definition 5.1. An intersecting pair with conic points is a subset, W, of 
sc T*{X)\qx which is a union of the closure of a smooth Legendrian subman- 
ifold, W , and which is a finite union of global sections of the form 

W^(Xj) = {(y, 0, Xj)}, and contains W\W. We also require W to have an 
at most conic singularity at fx = 0; that is, it is smooth if polar coordinates 
are introduced along W \ W. 

The process of introducing polar coordinates along W \ W can be given 
an invariant meaning and is then called blow-up. We denote the blown-up 
manifold by W. 

The metric g induces a metric h on the boundary as nearby it is of the 
form 

t 2 + h'(y,n) +xg 
as a function on SC T*X, we obtain h! from h via the isomorphism 

dy , 
/i. — i — ► fi.dy. 

x 
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Example 5.1. For each y' G dX and ^ A G R, let Gy(A) be equal to the 
set of (t, y, //), such that r 2 + |/i| 2 = A 2 ,/x 7^ 0, and putting /i = such 
that 

r = |A| cos(s) 

|A*I = |A|sin(s) (5.1) 
= exp(sHi h )(y',fi') 

where s G (0,vr), G T*3X, and £') = 1. Then G^(A)U {(A, y, 0)} 

is an intersecting pair with conic points. We denote this pair G{\). This is 
the pair in which we are interested here. The set G"(A) = {(— A, y, 0, y' , 0)} 
is also important in our construction. Note that G"(A) is the initial or 
incoming surface and that Gr(— A) is the outgoing surface. Note that in the 
coordinates defined by (5.1), G"(A) is s = 7r and A) is s = 0, when A is 
positive. 

Associated with these intersecting pairs at each conic point is a unique 
homogeneous Lagrangian submanifold A(W,Xi) of T*{dX). For the pair 
G(X) we are interested in, this is precisely the relation of being n apart 
along a lifted geodesic. (See Proposition 4 of |]15[ .) For simplicity, we shall 
henceforth take A to be positive. The A negative case is similar, or could be 
deduced from the positive case. 

Melrose and Zworski associated to any such intersecting pair a class of 
smooth functions whose asymptotics on approach to the boundary are de- 
termined by symbols on the Legendrians. A symbol bundle over the smooth 
Legendrian W(A) in the pair W can be defined and is denoted E m ' p . The 
sections of this bundle are of the form 

aS p - m \dx\ m - n ^ 
1 

where a is a smooth section of C°°(W; £l£ <8>M^), S is a defining function of 

1 

the boundary of W, M is the Maslov bundle, and is the b— half density 
bundle. For G above, one could take S = sins. Melrose and Zworski remove 
this singularity at the endpoints by rescaling but for us it will be easier not 
to do so. 

Proposition 5.1. If W{\) is an intersecting pair with conic points then 
there is a class of smooth half- densities on X° , denoted I™' P (X, W), such 
that P| /™' P (X, W) is equal to the class of half- densities vanishing to infinite 

m,p 

order at the boundary. There exists a symbol map 

v sc , m , P : I™' p (X,W, sc n 1/2 ) - C°°(TM m ' p ) 
which gives a short exact sequence 

-» I™ +1 ' p (X,W, sc n 1/2 ) -» ITc P {X, W, SC Q 1/2 ) -> C°°(W;E m ' p ) -» 0. 
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This is Proposition 12 from [El|. The Legendrian half-densities of order 
m are given locally away from the conic points by oscillatory integrals 

u=(27r)-?-5 [ e i ^ y ^a(x,y,u)x m - k / 2+n / 4 du, (5.2) 



with a smooth on [0, e) x U x U' with U, U' open and (j> parameterizes W, 
that is 

W = {(y,-^(y,«),d tf 0):d u ^ = O}. (5.3) 

Near the conic singularity a more general form is required and we refer 
the reader to [15|. The order m here is adjusted by — k/2 + n/4 and so 
is inconsistent with section |3|. We keep this inconsistency as the second 
definition provides for better invariance properties but less clarity. 

An important related fact we need to know is how Legendrian distribu- 
tions map under scattering pseudo-differential operators. We recall Propo- 
sition 13 from [15|. 

Proposition 5.2. Suppose P € ^ l j k (X, ^fi 1 / 2 ) has symbol x k pk + x k+1 pk+i 
with respect to a product decomposition of X near dX, and suppose that 

wc sc n x (x) 

is a smooth Legendre submanifold. Then for any m 6 R, 

P : ITc(X, W; SC Q 1/2 ) -» I™ +h (X, W; sc 1/2 ) (5.4) 
a sc>m+k (Pu) = (pk\G) a sc,m(u) ® \dx\ k . (5.5) 
Furthermore, if Pk vanishes identically on W then 

P : T£{X, W; ^n 1 / 2 ) -» I™ +k+1 (X, W; sc VL l > 2 ) 

and 

o- S c,m+k+\{Pu) = 
1 /. , n 



L V + ( I(* + 1) + m - ^) + p k+1]w ) a ® |dxr +fc+1 -f 



where c SCjm (u) = a<8> |dx| m * and V is i/ie rescaled Hamiltonian vector field 
associated to p k . 



We omit the definition of the rescaled Hamiltonian vector field but recall 
that for A on the pair G we are studying it is equal to 

2 A sin s— 

OS 

in the semi-global coordinates given by (5.1). 

We also need two push-forward theorems, Propositions 16 and 17, from 
[fl5| . They relate the singularities of the scattering matrix to the asymptotics 
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in small x of the Poisson operator. Given a product decomposition near the 
boundary, there is a natural pairing 

B -. c-™{x, sc q}-/ 2 ) x c^idx-^n 1 / 2 ) -> c~°°{[o,e), sc n 1/2 ) 

(5.6) 

n— 1 f 

B(u,f) = x— / u(x,y)f(y). (5.7) 

J dX 

Proposition 5.3. For any intersecting pair of Legendre submanifolds with 
conic points, W, the partial pairing (5.7) gives a map 

B : I™> P (X, W; SC Q 1/2 ) xC°°(dX; sc ^/ 2 ) ^ £ I p+ ^ ([0, e ), W' (ff^Q 1 / 2 )) 

j 

where the W'(fj) = {(0, — Tjdx/x 2 )} are the Legendre submanifolds corre- 
sponding to the components of W# and 



B(u,f) = J2e- if/x x p+n / 4 Q° f .(u,f) 

with 



3 

j 



dx 



+ 0{x 



p+n/4+l\ 



QUu)er ph ; 4 (dx,A(w,r)), 

and the principal symbol of Q%{u) determines and is determined by the lead 
singularity of the principal symbol of u on W on approach to W'(fj). 

When the Legendrian distribution is actually associated to a smooth Leg- 
endrian submanifold the push-forward becomes much simpler and this sim- 
plifies the construction of the Poisson operator. 

Proposition 5.4. If G is a smooth Legendre variety and u E I™(X, G 1 , sc VL l l 2 ) 
near r = r, then the distribution Q% is a Dirac delta distribution. 



6. The Poisson operator in the general case 

In this section, we apply the calculus reviewed in Section || to construct 
the Poisson operator and prove that the higher order scattering matrix is 
indeed a zeroth order, classical, Fourier integral operator, proving Theorem 



1.2 and Corollary 1.1 . We shall refer heavily to [15], Section 15 as our 



construction is a modification of the one there. 

We assume a product decomposition of X close to the boundary of the 
form (|1.1| ) has been chosen and is fixed throughout this section. We then 
have as in |l5| that A, the intrinsic Laplacian acting on scattering half- 
densities on X, induces an operator 

A x G Diff 2 c (X x dX, sc VL l l 2 (X x dX)) 
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by 



A; 



dx 
x 2 



dy 

r.n—1 



dy' 

■y,n— 1 



A «(.,!/) 



X 2 



dy 

r.n—1 



dy' 

T,n— 1 



where (x, y, y') is a point in XxdX. Throughout this section V will be a short 
range higher order perturbation, that is a scattering differential operator of 
order 2k — 1 as a differential operator and order 2 at the boundary. 

We recall from , using the notation of that 

Lemma 6.1. The symbol at the boundary of Ax is p = Po + xp\ with po = 
t 2 + h(0, y, n) and with p\ equal to —i(n — l)r + c where c = -§^h(x, y, [i)\ x =o 
is quadratic in \i. 

Our choice of product decomposition ensures that there is no t term in 
c. We need to compute the symbols at the boundary of — X 2k + V. The 
short range assumption on V ensures that it does not contribute - indeed 
this is one reason why this definition of short range is appropriate. Now we 
can decompose 



(A fc - X 2k ) = Q(X)(Ax - A 2 ), 



with 



fc-i 



Q(A) 



.7=0 



fc-1- 



We deduce that the lead symbol of (A^ — \ 2k ) is Pq — \ 2k and the second 
symbol (subprincipal term) is equal to k\ 2 ^ k ^^ times p\ on the zero set of 
the principal symbol. It also follows that the Hamiltonian on the boundary 
of the lead term is just A;A 2 ^' -1 ^ times that of po on the zero set of the 
principal symbol. The fact that both these terms have been changed simply 
by multiplication by kX 2 ^ k ~^ ensures the simplicity of extending results from 
the k = 1 case to the higher order case. 

We also note the following lemma which is important in our construction 
to show that the transport equations are solvable. 

Lemma 6.2. If L G I m '-^(X x dX, (5(A), sc r2 1 / 2 ) is such that 



then 



(A| - X 2k + V)L e I m+j 'i{X x dX,G{X), sc Q 1/2 ) 
(A k x - X 2k + V)L G I m+i 'i{X x dX,G{X), sc Q 1/2 ) 



This is a modification of Lemma 15 from 1£] and in fact, the k = 1 case, 
though not explicitly stated, is essential to the construction there also. 



Proof. The proof is no different from that of the special case in M n (Lemma 

H). □ 
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Proposition 6.1. For any / A g i there exists 

K G I m ' p (X x dX,G{\); sc ^ 1/2 ) 

such that 

(A k x - X 2k + V)K G C°°(X x dX- ^n 1 / 2 ) and 
Q°x(K) = Id, 

with 

2n-l 1 
m = —'P=—^ 

and the principal symbol of K on G is 

(sin s) 2 

where C(y,£i) is a non-zero smooth function. 

Proof. As in |§, we first construct K b G I m ' p (X x <9X, G(A); ^Q 1 / 2 ) such 
that 

(A* - \ 2k + V)K h G /J C (X x OX, G tt (-A)) and (6.1) 
Q°x(K b )= Id. (6.2) 
We construct if 6 as an asymptotic sum of 

2n-l | ■ _ 1 

i^ GJ sc 4 +J ' 3 (X x dX,G(X); sc n 1 ^ 2 ). 

We wish 

2n-l , r, 3 

(A^ — A + V)Kq G J sc 4 3 (X x aX,G(A); sc n 1/2 ) and 

(6.3) 

MQx(Ko)) = ^o(Id) (6.4) 
and then it is automatically in 

271-1 i r, 7 

I sc 4 + ' 3 (X x OX, G(A); SC VL 1 / 2 ) 



by Lemma p.2[ We also want 



(1 } \ 2n-l i ■ , n 3 

X^J 4 3 ^(Ixai,^);^ 2 ) 



and this, of course, implies that it will also be an element of 

2n— 1 i „■ i o 7 

J sc 4 +3+ ,I (IxM,G(A); sc H 1/2 ). 



Now near GnG'(A), where G is smooth, we can as in give an explicit 
construction and it is then only necessary to have that the principal symbol 
of Q x (Kq) is equal to 1 to ensure that Q x (K b ) = Id . We look for Kj of the 
form 

xjjWMyxa^y^y^X^^aj G C°°{X x dX), 
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with v a fixed scattering half-density and 4> the cosine of the Riemannian 
distance from y toy' . Let cij be the restriction of dj to x = 0. Taking geodesic 

normal coordinates, y, about each y' the transport equations for a.- is of the 
form 

(y ■ d y + j)d d + ftj-Oj- = cj G G°°(X x 8X) 

near y = where Co is identically zero and 6j vanishes quadratically at y = 0. 
So as in [15], the terms Kj exist sc-microlocally close to G"(A). 

We now need to continue each up to G*{— A). We do so by solving 
transport equations for the principal symbols and iteratively solving away 
the error. 

The principal symbol of Kq, a m (Ko), is of the form 

li 



\ds\z\dy\2 \dfi\- 



(sin s) 



■2 



On the lifted geodesic (3{s) the sub-principal term c(/3(s)) = 2k\ sm(s)d(f3(s)) 
for some smooth d. From Proposition 5.2, the transport equation for b is 

2£;A 2 ^ 1 / . , , d n-n\ , . . . ,,\ , 
f sin(s)- + f — ^— J cos(s) + zsm(s)(i(/3(s))J 6 = 

Writing b = ( e i f d Ws'))ds' gin^) 1 ^)^ we thus have 

— - 

ds 

This means that 

6 = Csin( S )V e -/^'))^'. (6 .5) 



n-l 
2 



As s — > 7T— , that is near G tt (A),this has a singularity of the form (vr — s) 
and as s — > 0+, of the form s 2 . 

As the order on G"(±A) is — i, the order p — m is equal to the order of 
singularity and thus the solution of the transport equation is a legitimate 
symbol and we can construct Kq. 



Now by Lemma |6.2j , 

(A*. + V - X 2k )(K ) G /-^+ 2 '-i+ 2 (X x 9X, G; sc fW 2 ), 
and we look for 

such that 

(a| + v — A 2fc )(^ + k x ) g j-^+ 3 >-*+\x x ax, G; sc ft 1/2 ), 
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±(e i S d ^ ds ' S m(s) 1 ^ + %)=g(t 



and thus by Lemma |0| is in /- 22 t j -+ 3 ~i+ 2 (X x dX, G; ^VL 1 / 2 ). Letting the 

2n — 1 I -i 2n — 1 

principal symbol of K\ be Oi|ax| 2 h * times the trivializing density 
above, we obtain a transport equation; arguing as above, it becomes 

as V / 

= g(s)e- l I d ^ ds 'sm(s) ]1 i 1 , 

with g(s) a smooth function on [0, tt] (and smoothly depending on the sup- 
pressed parameters). Canceling, we obtain that 

d_ 

ds 

which has a solution in the appropriate symbol class. The same argument, 
after appropriately shifting indices, constructs all the terms Kj. 

Asymptotically summing, we obtain K b such that 

(A k x + V- \ 2k )K h e / 7/4 (G«(-A)). 

These errors can now be removed by an iterative construction of their Taylor 
series, cf Lemma 16 of [fig] , and we obtain K as desired. □ 

We have therefore constructed the Poisson operator modulo smooth terms 
vanishing to infinite order at the boundary. We wish to remove this error 
by applying the resolvent. As before, this is possible even if there is embed- 
ded discrete spectrum, as the error is orthogonal to the eigenspace at that 
energy. This is seen by a simple integration by parts, since the elements 
of the eigenspace are Schwartz they are orthogonal to the image of smooth 
functions of tempered growth. We have thus constructed P(A), the Poisson 
operator for the problem, as a paired Legendrian distribution. The remain- 
der of the proof that S{X) is a Fourier integral operator now follows as in 



[|15| , Proposition 19. 



To prove Corollary |l.l| , one observes that if -Pi (A) is the Poisson operator 
for A — A 2 , then 

(A fc + V - X 2k )Pi(X) = Q(A)(A - A 2 )Pi(A) + VPi(X) = VP^X). 

(6.6) 

This means that the Poisson operator for A k + V — X 2k can be constructed 
as a perturbation of -Pi(A) and as VP\{X) G / ~ +^-4+^ the scattering 
matrices will agree to order 1 — I, and the corollary follows. Note, however, 
that even if one fixes the perturbation, the scattering matrices will not agree 
to order more than l — l as the solutions of the transport equations will differ. 
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